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Abstract. - We consider the tunnel current between two ferromagnetic metals from a per- 
spective similar to the one used in superconductor/superconductor tunnel junctions. We use 
fundamental arguments to derive a Josephson-like tunnel spin current /-^p'" qc sin(&i — 62). 
Here the phases are associated with the planar contribution to the magnetization, (cjcj) ~ e'^. 
The crucial step in our analysis is the fact that the z-component of the spin is canonically 
conjugate to the phase of the planar contribution: [6, S'] = i. This is the counterpart to the 
commutation relation [ip, iV] = z in superconductors, where ip is the phase associated with 
the superconducting order parameter and A'^ is the Cooper pair immber operator. We briefly 
discuss the experimental consequences of our theoretical analysis. 



In recent years the possibility of using the spin degrees of freedom in addition to the 
charge degrees of freedom in electronic devices has open new perspectives in condensed matter 
physics. Indeed, the field now known as spintronics or magnetoelectronics is rapidly developing 
and there are even commercial devices available [1]. If the current density due to an electron 
of spin fj is jo-, then the usual electric current is just the sum of the electronic currents for 
each spin degree of freedom, 'f^'^'^s'^ = jo-. The spin current density, on the other hand, can 
be defined by j**?'" = crjo-- Although such a definition of spin current density is used very 
often, it should be kept in mind that that the true spin current density is actually a tensor 
and not simply a vector like the charge current density. Indeed, if M is the magnetization, 
the continuity equation for the spin current is given by dMi/dt -f df^J^'^™ = 0, where J^"^'" is 
the spin current density tensor. 

The most familiar situation in the theory of electronic transport corresponds to the case 
where a charge current flows in the system and the spin current is zero. However, the opposite 
situation is also possible in some ferromagnetic systems, i.e., a spin current flowing in the 
absence of any charge current [2-4] . A relevant problem to study in this context is the behavior 
of the spin current across a tunnel junction between two ferromagnetic metals [5], especially 
whether there is phase coherence of the spin current across the junction. Phase coherent 
charge transport is a well known feature in tunnel junctions between superconductors and is 
responsible for the Josephson effect [6]. In this paper we will study a similar effect for a spin 
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current in ferromagnet/ferromagnet (FM/FM) tunnel junctions. Recently it was pointed out 
that a charge Josephson-like effect would be possible in a spin valve [7] , but our analysis below 
does not confirm such a prediction. The spin Josephson effect in itself is not new and was 
observed some time ago in '^He-B weak links [8,9]. The observation of this effect establishes the 
existence of spin supercurrents in '^He-B. It is plausible that a similar effect may occur in other 
systems, like in tunnel junctions involving ferromagnetic metals. Another possibility is to use 
thin film helimagnets, where spin transport without dissipation was discussed recently [10]. 

The main result of this letter will be a microscopic derivation of the spin Josephson effect in 
FM/FM tunnel junctions. The key observation in our analysis is that in a ferromagnetic metal 
spin flip averages like (c|c|) are nonzero and play a role analogous to the superconducting 

order parameter (cjcj). This same analogy was used recently in discussions of a Josephson- 
like effect [11-14] observed in double layer quantum Hall pseudo-ferromagnet [15]. There the 
spin label is associated with the layer such that an analogy with exciton superfluids [16] is 
also possible [17]. 

In order to make our analysis of the spin Josephson effect self-contained, we first recall some 
basic facts about superconductor/superconductor (SC/SC) tunnel junctions. There the phase 
plays a crucial role leading to the celebrated Josephson effect [18]. Indeed, the Josephson 
effect is a consequence of the fact that the phase (f of the order parameter is canonically 
conjugate to the Cooper pair number operator N, [ip,N] = i, which implies the uncertainty 
relation AipAN ~ 1 [18]. The superconducting current is given by Isc = 2e((7V2 — Ni)), 
where e is the electric charge and iVi and N2 are the Cooper pair number operators for the 
superconductors in each side of the tunnel junction. Using then the Hamiltonian for the tunnel 
junction 

H = -Ej C0S(^1 - ip2) + ^(^1 - iV2)2, (1) 

where C is the capacitance of the tunnel junction, and treating N.^ and ipi as classical 
canonically conjugate variables, we easily obtain I^j'^ = {Ej /e) sin Aip, Aip — 2eV, where 
A(p = (pi - (P2 and F = 2e(iVi - iV2)/C [19]. 

In the case of FM/FM tunnel junctions it is the z-component of the spin that is canonically 
conjugated to the phase: [9, S^] = i. This follows from writing the raising and lowering spin 
operators as [20] 




5"+ = exp(i6 

S- = ^ + § j '[^^^2) '^'^P(-*^)- (2) 

Thus, we obtain the uncertainty relation A^A^*^ ~ 1. Furthermore, we can derive the commu- 
tation relation [S^ , e='=*'"^] = ±e^""^ which in turn implies the commutation relations defining 
the algebra of spin operators. Using the representation Q we can establish a classical Hamil- 
tonian for the FM/FM tunnel junction out of which a spin current /-^p'" — fisiSf^ — Sf) is 
derived: 



H 



-EjS' 



^^^^ 



SI 



(3) 
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Table I - Analogies between SC/SC and FM/FM tunnel junctions 



SC/SC junction FM/FM junction 

Voltage V Spin voltage 14 

/sc = (£;j/e)sinAv9 r/" = {2E.jS^ / ^iB)s■^n M 



where for convenience we have introduced the "spin capacitance" C^, which will be specified 
later in the discussion of the microscopic origin of the effect. The coeffcient EjS^ of the first 
term of Eq. Q corresponds in fact to an exchange energy between the magnetizations in 
both sides of the tunnel junction. The value of EjS^ will be determined later in terms of 
parameters of the microscopic model. 

Using the classical Hamiltonian equations of motion we easily derive: 

M = 2fiBVs, if"' - {2EjSyfiB) sinAe, (4) 

where we have defined the spin voltage as Vg — {^b / Cs){Sf^ — S^). Thus, a spin current in 
a FM/FM tunnel junction would behave in the same way as the superconducting current in 
the SC/SC tunnel junction, see Table L 

Now that we have motivated the problem on the basis of a straightforward physical analogy, 
we are ready to derive microscopically the spin Josephson effect in FM/FM tunnel junctions. 
Before starting, it is worth to compare once more the role of the phases in superconductors 
and ferromagnets. In superconductors a phase transformation of the electron destruction 
operators involves the particle number operator: 



i — (n 



exp 
e*'^/2 c,(0). 



Ccr(O) exp 



■"fit 



1) 



(5) 



where n — clyCa and from the first to the second line we have just solved the equation of 
motion dca/d^p = i[{n — l)/2,Co-]. Similarly, for a ferromagnetic metal the relevant transfor- 
mation is 



exp 



i — (nt 
2^ ^ 



Co.(0) exp 



2^ ^ 



CaiO). 



(6) 



Eq. ^ is thus associated to gauge transformations Ca = e^'^^^c^r out of which we obtain that 
(cjcj) = e~*'^(c|cj). Eq. ®, on the other hand, yields (cjq) = e~*^(c|c4). 

We will consider the ferromagnetic metals in each side of the junction in the simple frame- 
work of the Stoner model for itinerant ferromagnetism [21]. The interaction is given simply 

by 



Hr 



U / d^r n|(r)n4(r). 



(7) 



which can be rewritten as 
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H'j = Jd'rl^^[m\r) + \A{r)\']~c\{r)c,{r)A{v) 

~ At(r)c|(r)cT(r)-m(r)K(r)-nx(r)]}, (8) 

upon a Hubbard-Stratonovich transformation. It is easy to see from the equations of motion 
for the auxihary fields that mean-field theory corresponds to 

At = (2C//3)(5+), A = (2C//3)(5_), m = (2C//3)(5,), (9) 

where the spin operators are defined in terms of fermion operators as 5*+ = c|c| and 5*2 = 
K -'^i)/2- 

The electron field operators in the left and right subsystems will be denoted as Ca and 
/ct, respectively. The bosonic auxiliary fields in the left and right subsystems will be distin- 
guished respectively by labels L and R and assumed to be uniform in a mean-field approxi- 
mation. Futhermore, for simplicity we will also assume that the left and right ferromagnetic 
metals are the same. Thus, we will take m| = = vr? and |Ai| = |Aj^| = |A|, but 
in general 7^ w/j, since we can have m^, = toq and = —mo. Also, the phases 
of Ai and A^ are generally not the same. The tunnel Hamiltonian is given as usual by 
Ht — X^k p X]cr(^kpCk£^/pCT+b.c.), where Tkp is the tunnel matrix element. The charge current 
is given by /^^argc ^ _e(iV^(t)), where iV^ = Ek.. cL^k. and iV^ = i ^ E.(^^kpcL/p<x - 
^kp/paCka), while the spin current is given by /"P" = -hb{SI) where S'£ = Y.k,a '^4^aCko- 

and SI = i J2w,p J2a (^iT^pC^fpa ~ T*p/t^Ck<,). 

We will evaluate the spin current using linear response theory. In this approximation we 
have /spi" = 2/isIin[S(ey -I- iS)], where 6 0+. In terms of the Matsubara formalism we 
have 



1^1 k,p 



<jG1{k, iuji)G'^{p, iuji — ioj) 

(7 

The Green functions appearing in Eq. (|10|l are given in the mean-field approximation by 

nan ■ s iuj-ei^ + cmiL 

= ^ (11) 

[luj — ek) — TO^ — |A|"^ 

with similar expressions for the right side Green functions. The dispersion energies are ek = 
Sk — ^j^■L and — — fij^, where fiL — /^j?, = eV. The spin current can be separated into two 
parts, one due to direct spin current tunneling and another one due to spin flip tunneling. 
It is the latter that is responsible for the spin Josephson effect and arises from the product 



(iuj — 


ek)^ — ra^ — 


|AP 








{iuj — 


ek)^ — — 


|AP 




Al 




{iuj — 


ek)^ — rn^ — 


|AP 
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between the F'j^^j^ Green functions. The direct spin current vanishes at zero voltage while the 
spin flip current does not. Explicit evaluation gives 

where Ad is the phase difference between the left and right subsystem as before. The function 
S{a, b) is given by 

f{v + h) + f{v -b)- f{u + b)- f{u - b) 



S{a, b) = j du j dv p{u) p{v) 

fin, -h\ - fill -I- Ml 

(13) 



a + V — u 

f{v + b)-f{u-b) f{v-b)-fiu + by 



a + 2b + V — u a — 2b + V — u 

where p(e) is the density of states and /(e) is the Fermi-Dirac distribution function. The 
function S{a,b) satisfy the properties S{a,0) = and S{0,b) / 0. Eq. ifT^ establishes 
microscopically the spin Josephson effect in FM/FM tunnel junctions. The Josephson-like 
spin tunnel current is nonzero only if |A[ ^ 0. Furthermore, it is maximum for m = 0, which 
corresponds to the case of a tunnel junction between planar ferromagnets. 

From Eq. (|12|l we can easily determine the exchange magnetic energy across the junction, 
since ly^™ = dE/dA9, where E is the energy of the tunnel junction. In this way we obtain 
that the coefficient EjS^ appearing in Eq. Q is given by the coefficient of the sine of Eq. 
(O with F = 0. 

The spin Josephson effect as given by Eq. (jl2|l implies a spin current even in the absence 
of voltage. Now we will show that when tul — rnji is nonzero the above effect is actually an 
AG-like effect. In order to show this we have to use a gauge transformation based on Eq. 
(O in the Lagrangians for the left and right subsystems, i.e., we set c^a = e^^^^^^Cko- and 
/ko- — e^'^^"^^fka- For instance, we obtain for the left subsystem 

Lr = E U - ld,e, + am, + ^) 5^. - fff^ (14) 

with a similar expression for the right subsystem. In analogy with the Josephson effect in 
superconductors, we can compensate for the phase transformation by writing dfO, = 2mi, 
and OiOr = 2mfl, out of which we obtain 

= 2{mL - mil) = ^fisV,, (15) 

where we have defined a "spin voltage" Vs = (m^ — mji)/fiB- Using the expression for m in 
Eq. for the left and right systems and comparing Eq. IjlSI with Eq. 10}, we obtain that 
the "spin capacitance" is given by 

ft. = M. (16) 

It is worth to discuss once more the similarities between FM/FM and SG/SG tunnel junc- 
tions. In the latter the time evolution of the phase difference is proportional to the average 
particle number difference [22], while in the former it is proportional to the difference of 
the ^-projection of the magnetization across the junction. Note that in superconductors a 
Hubbard-Stratonovich transformation similar to the one in Eq. (|SJ| holds with the spin oper- 
ators replaced by the pseudo-spin operators 5+ — c|c|, S- — 5^, = {n^ + — l)/2. 
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If m| = mfj. = 777?, the AC-like spin Josephson effect occurs only if the left and right 
subsystems have the spin projection on the z-axis anti-parallel to one another. By solving Eq. 
H15|l with rriL = —mn — m and using Eq. (|12|l we obtain the important result 

= 2p'+|S) ^^'^' v/^^^^T^) sin(A0o + irnt). (17) 

The spin Josephson effect can be probed by trying to detect the induced electric fields 
associated with it [23]. It can be seen from Maxwell equations that a spin current should 
induce an electric field in the same way that a charge current induces a magnetic field. Indeed, 
in absence of voltage the induced electric field should satisfy 

d^Ey - dyE., = ~A^^^lB^{S,)/^t - Airlf'', (18) 

where we have assumed that dH/dt = 0. 

In conclusion, we have shown that a spin Josephson effect should occur in FM/FM tunnel 
juntions. In the more general case where the FM order parameter is three-dimensional, the 
effect is necessarily an AC-like one, with the oscillations in time being associated to the z- 
component of the magnetization. Finally, we would like to stress that our analysis is completely 
different from Josephson effects discussed previously in FM/FM junctions where a charge 
current is discussed [7]. Here we have considered the phase coherence of a spin tunnel current. 
Our point of view is similar to the one of Ref. [10] where averages like (cjcj^) also play a 
crucial role. In Ref. [10] the phase coherence arises from the fact that a thin film helimagnet 
is being considered while in our case it is due to the existence of a boundary between two 
ferromagnetic metals. 

* * * 
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